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We present a global fit for a s (mz), analyzing the available C-parameter data measured at 
center-of-mass energies between Q = 35 and 207 GeV. The experimental data is compared to a 
N 3 LL' + 0(a 3 ) + Di theoretical prediction (up to the missing four-loop cusp anomalous dimen¬ 
sion), which includes power corrections coming from a field theoretical nonperturbative soft func¬ 
tion. The dominant hadronic parameter is its first moment f2i, which is defined in a scheme which 
eliminates the O(Aqcd) renormalon ambiguity. The resummation region plays a dominant role 
in the C-parameter spectrum, and in this region a fit for a s (mz) and Hi is sufficient. We find 
a 3 {mz) = 0.1123 ± 0.0015 and Di = 0.421 ± 0.063 GeV with x 2 /dof = 0.988 for 404 bins of data. 
These results agree with the prediction of universality for Hi between thrust and C-parameter within 

1-(T. 


I. INTRODUCTION 


In order to study Quantum Chromodynamics (QCD) 
accurately in the high-energy regime, it is useful to ex¬ 
ploit the wealth of data from previous e + e~ colliders 
such as LEP. Here the final states coming from the under¬ 
lying partons created in the collisions appear as boosted 
and collimated groups of hadrons known as jets. Event 
shapes have proven to be very successful to study these 
collisions quantitatively. They combine the energy and 
momenta of all of the measured hadrons into an infrared- 
and collinear-safe parameter which describes the geomet¬ 
ric properties of the whole event by a single variable dis¬ 
tribution. Due to their global nature event shapes have 
nice theoretical properties, making it possible to obtain 
very accurate theoretical predictions using QCD. Most 
e + e _ event shape variables quantify how well the event 
resembles the situation of two narrow back-to-back jets, 
called dijets, by vanishing in this limit. Because the dijet 
limit involves restrictions that only allow collinear and 
soft degrees of freedom for the final-state radiation, such 
QCD predictions involve a number of theoretical aspects 
that go beyond the calculation of higher-order pertur¬ 
bative loop corrections. These include factorization, to 
systematically account for perturbative and nonpertur¬ 
bative contributions, and the resummation of large log¬ 
arithmic corrections by renormalization group evolution. 
Comparisons of predictions for event shapes with experi¬ 
mental data thus provide non-trivial tests of the dynam¬ 
ics of QCD. 

Due to the high sensitivity of event shapes to jets 
induced by gluon radiation they are an excellent tool 
to measure the strong coupling a s . For more inclusive 
hadronic cross sections (like e + e _ —► hadrons) the a s 
dependence is subleading because it only occurs in cor¬ 
rections to a leading order term, while for event shapes 
the a s dependence is a leading-order effect. For this rea¬ 
son, the study of event shapes for determining a s has 
a long history in the literature (see the review [Tj and 
the workshop proceedings ®), including recent analyses 


which include higher-order resummation and corrections 
up to 0(ctg) [3VF12]. 

Several previous high-precision studies which deter¬ 
mine a s (mz ) [HI51I5HTT] focus on the event shape called 
thrust T3j, 


r = 1 — T = min 1 — 


E j 

Ej \Pj\ J ’ 


(1) 


where n is called the thrust axis and it follows from the 
above equation that 0 < r < 1/2. Another event shape, 
known as C-parameter El US], can be written as: 

c = 3 E i,j IPillPjlsin 2 % 

2 (Ei Wi\f 

where Oij gives the angle between particles i and j. It 
is straightforward to show that 0 < C < 1. In a pre¬ 
vious paper ]12] we computed the C-parameter distribu¬ 
tion with a resummation of large logarithms at N 3 LL' 
accuracy, including fixed-order terms up to 0(a 3 ) and 
hadronization effects using a field-theoretic nonperturba¬ 
tive soft function. These results were achieved by using 
the Soft Collinear Effective Theory (SCET) [TBTE0] . Our 
results for C are valid in all three of the peak, tail, and 
far-tail regions of the distribution, and are the most ac¬ 
curate predictions available in the literature, having a 
perturbative uncertainty of ~ 3% at Q = mz for the re¬ 
gion relevant for a s (mz) and Hi fits. The same accuracy 
was previously achieved for thrust, where the remaining 
perturbative uncertainty in the r distribution is ~ 2% in 
this region [5], In this paper we make use of these new 
C-parameter theoretical results [T2] to carry out a global 
fit to all available data, comparing the results with the 
analogous global fit for thrust [9] where appropriate. 

Since both r and C vanish in the dijet limit, it is worth¬ 
while to contrast them in order to anticipate differences 
that will appear in the analysis. Differences between C 
and r include the following: 
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a) calculating r requires identifying the thrust axis 
with a minimization procedure, while C does not 
involve a minimization; 

b) r has a single sum over particles while C has a 
double sum; 

c) the size of the nonperturbative region, where the 
entire shape function is important, is larger for C 
compared to r due to an enhancement by a factor 
of 37 t/2; 

d) the resummation region for C is larger than that 
for thrust since the logarithms appear as ln(C/6) 
compared to ln(r), which increases the range of C 
values that are useful for a s fits; 

e) fixed-order predictions for the thrust cross sec¬ 
tion are smooth across the threshold where non- 
planar events first contribute, r = 1/3, while the 
fixed-order C-parameter cross section has an inte- 
grable singularity at this threshold, C/hcmMer = 3/4. 
The singularity for C comes from the fact that 
the leading-order distribution is not continuous at 
C = Cghouider- (The C'-distribution can be made 
smooth here using LL resummation |2lj.) 

A key similarity is that in the dijet limit ( C , r <C 1) the 
partonic cross sections for thrust and C-parameter are 
related up to NLL by using tnll = Cnll/ 6 US- This 
relation quantifies several qualitative relations between 
C and t in the dijet region. 

Recent higher-order event-shape analyses 01 IMB 
have found values of a s (mz) significantly lower than the 
world average of a s (mz) = 0.1185 ± 0.0006 [23] which 
is dominated by the lattice QCD determination pH] . 
This includes the determination carried out for thrust 
at N 3 * LL' + O(a^) in Ref. [5]R which is also consistent 
with analyses at1 N 2 LL + O(o“ ) [HI El] which consider 
the resummation of logs at one lower order. In Ref. [25] 
a framework for a numerical code with N 2 LL precision 
for many e + e - event shapes was found, which could also 
be utilized for a s fits in the future. In Ref. <9J it was 
pointed out that including a proper fit to power correc¬ 
tions for thrust causes a significant negative shift to the 
value obtained for a s (niz), and this was also confirmed 
by subsequent analyses m Recent results for a s (mz) 
from r decays mi, DIS data [283, the static potential for 
quarks [25] , as well as global PDF fits [351 El] also find 
values below the world average. With the new analysis 


1 Note that results at N‘ 5 LL require the currently unknown QCD 

four-loop cusp anomalous, but conservative estimates show that 

this has a negligible impact on the perturbative uncertainties. 

Results at N 3 LL' also technically require the unknown 3-loop 

non-logarithmic constants for the jet and soft functions which 
are also varied when determining our uncertainties, but these 
parameters turn out to only impact the peak region which is 
outside the range of our ct s ( m z) fits in the resummation region. 


we present here, we provide another event-shape determi¬ 
nation of a s (mz ) with a high level of precision. We will 
also simultaneously examine the leading power correction 
to the distribution, which should be universal between 
thrust and C-parameter. 

This paper is organized as follows. In Sec. [IT] we review 
the theoretical calculation of the C-parameter cross sec¬ 
tion, presented in more detail in Ref. [12]. The details on 
the experimental data and fit procedure used in our anal¬ 
ysis are given in Secs. [ITT] and [IV] In Sec. [V] we present the 
results of our fit for a s (mz) and the first moment of the 
nonperturbative soft function, O/'. Fits which include 
hadron-mass effects are discussed in Sec. IV El In Sec. IVT1 
we make predictions for the peak and far-tail regions of 
the distribution, which are not used in our fit, and com¬ 
pare those regions to experimental data. The universality 
of fii is discussed in Sec |VII[ where we compare our re¬ 
sults with the previous fit done using thrust in Ref. [5]. 
Finally, Sec. VIII contains our conclusions. We also in¬ 
clude three appendices. The first, Appendix [A] contains 
the formulae needed to calculate the profile functions, 
the second Appendix [B] contains results that support our 
choice for the definition of the renormalon free pa¬ 
rameter to use for C ; and the third compares fit results 
for thrust with our profiles and those from Ref. (5J. 


II. THEORY REVIEW 

Until a few years ago, the theoretical uncertainties re¬ 
lated to perturbative contributions as well as hadronic 
power corrections were still larger than the experimental 
uncertainties. The situation on the theory side has dra¬ 
matically changed with the calculation of O(a^) correc¬ 
tions m [71 [32HM] , and the pioneering use of SCET to ob¬ 
tain higher-order perturbative corrections in Ref. 0] , and 
to obtain accurate predictions for the full spectrum and 
a precision description of power corrections in Ref. [9]. 
In this section we review the theoretical work behind our 
calculation of dtr/dC, presented in Ref. [15] , 

SCET separates the physics occurring at the different 
energy and momentum scales relevant to the underlying 
jets whose properties are characterized by an event shape. 
For the C-parameter distribution the relevant scales are: 
(i) the hard scale /i#, which is related to the production 
of partons at short distances and is of the order of the 
center-of-mass energy Q ; (ii) the jet scale [ij ~ Q^/C/d>, 
which governs the formation and evolution of the jets; 
and (iii) the soft scale ~ QC / 6, which is the scale of 
large-angle soft radiation. All three scales are widely sep¬ 
arated in the dijet region C<1, where most of the events 
occur and where the distribution is maximal. There is 
one function associated to each one of these scales in 
the factorization theorem that describes the dominant 
contribution of the C-parameter distribution in the dijet 
limit: (i) the hard function H (the modulus squared of 
the QCD-to-SCET matching coefficient), which is com¬ 
mon to all dijet event-shape factorization theorems; (ii) 
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the jet function J T (given by matrix elements of quark 
fields with collinear Wilson lines), which is common for 
C-parameter [T2], thrust [5] and Heavy-Jet-Mass (p) 0; 
and (iii) the soft function Sg (defined by a vacuum ma¬ 
trix element of purely soft Wilson lines), which in general 
depends on the specific form of the event shape. Whereas 
the former two are perturbative (ph,Pj Aqcd)j per¬ 
mitting the calculation of the hard and jet functions as 
an expansion in powers of a s , the soft function has per¬ 
turbative corrections (ps ^ Aqcd) as well as nonper- 
turbative contributions that need to be accounted for 
(ps // Aqcd)- Renormalization evolution between the 
three scales pn, Pj , and ps sums up large logarithms 
to all orders in perturbation theory. It turns out that 
the soft function anomalous dimension for C and r are 
identical m, providing a connection between these two 
event shapes at every order of perturbation theory. Only 
corrections related to non-logarithmic terms in their soft 
functions, and the associated towers of logarithms, differ 
between C and r. 

The soft function can be further factorized into a par- 
tonic soft function 5g, calculable in perturbation theory, 
and a nonperturbative shape function Fc, which has to 
be obtained from fits to data. In the strict MS scheme 
this factorization was achieved in Refs. [551157], (Ana¬ 
lytic power corrections for the C-parameter distributions 
have also been studied in other schemes and frameworks, 
see e.g. Refs. [551440] . 1 

The treatment of hadronic power corrections greatly 
simplifies in the tail of the distribution, defined by 
QC 37tAqcd> where the shape function can be ex¬ 
panded in an OPE. Here the leading power correction is 
parametrized by Of, the first moment of the shape func¬ 
tion. The main effect of this leading power correction 
is a shift of the cross section, da(C) —> d a(C — Of/Q). 
Interestingly, this matrix element is related to the corre¬ 
sponding one in thrust by 


This relation was first derived using dispersive models 
with a single soft-gluon approximation in Ref. [JD] ■ The 
equality can actually be derived to all orders in QCD just 
using quantum field theory {41], but ignoring hadron- 
mass effects [55]. These hadron-mass effects can also 
be formulated purely with quantum field theory oper¬ 
ators [43] . Although they may in general give large cor¬ 
rections, hadron-mass effects turn out to violate Eq. ^ 
at only the 2% level, which is well below the 15% level 
determination of Hi that we will achieve here. When pre¬ 
senting the results of our fits, we parametrize the power 
correction using Hi defined in Eq. 0 to ease compari¬ 
son with our previous analysis which determined this Hi 
based on thrust mm- 


A. Factorized Cross Section Formula 

In order to understand the perturbative components 
of the C-parameter cross section we make use of the 
C-parameter factorization formula. To make the connec¬ 
tion to thrust simpler we will often use functions defined 
with the variable C = C/6. For the perturbative cross 
section we find [T2| : 

vM = f««3(4) 

Here J T is the thrust jet function, obtained by the con¬ 
volution of the two hemisphere jet functions, J T = J (g) J. 
J T describes the collinear radiation in the direction of the 
two jets. Its definition and expression up to O(a^) [4H - 
146] can be found in Refs. 13 El- The three-loop non- 
logarithmic coefficient of this jet function, j' 3 , is not 
known, and we vary it in our scans. The anomalous 
dimension of J T is known at three loops, and can be ob¬ 
tained from Ref. m- 

The hard factor H contains short-distance QCD effects 
and is obtained from the Wilson coefficient of the QCD- 
SCET matching of the vector and axial-vector currents. 
The hard function is the same for all event shapes for 
massless quarks, and its expression up to O(a^) [481454] , 
can be found in Ref. jj]. The full anomalous dimension 
of H is known at three loops, 0(a^) [451 f5Tl [55] . 

The soft function Sq describes wide-angle soft radia¬ 
tion between the two jets. It is defined as 

S d (l,p) = -l- { 0\trY T n Y n 6(e-^y^Y* n \0) , (5) 

where C is an operator whose eigenvalues on physical 
states correspond to the value of C-parameter for that 
state: C\X) = C(X) \X). Since the hard and jet func¬ 
tions are the same as for thrust, the anomalous dimen¬ 
sion of the soft function has to coincide with the anoma¬ 
lous dimension of the thrust soft function. Hence one 
only needs to determine the non-logarithmic terms of 
the C-parameter soft function. In Ref. m we computed 
it analytically at one loop, sf = — 7 r 2 CV/ 6 , and used 
EVENT2 to numerically determine the two-loop non- 
logarithmic coefficient s^, with the result 

sf =-43.2 ± 1.0. (6) 

The three-loop non-logarithmic coefficient of the 
C-parameter soft function, s^, is currently not known, 
and we estimate it with a Pade approximation, assigning 
a very conservative uncertainty. We vary this constant 
in our scan analysis. The precise definitions of j 3 and 
as well as can be found in Eqs. (A10) and (A12) of 
Ref. 121. 

In Eq. 0 the hard, jet and soft functions are evalu¬ 
ated at a common scale p. If fixed-order expressions are 
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used for these functions, then there is no scale choice that 
simultaneously minimizes the logarithms for these three 
functions. One can instead renormalization-group evolve 
from p to the respective scales pn ~ Q, I^J ~ < 51 / 0/6 
and ns ~ QC/ 6 at which the logs in each of H , J T , and 
Sq are minimized, and only use fixed-order expressions 
for these functions at these scales. In this way, large logs 
of ratios of the scales are summed up in the renormaliza¬ 
tion group evolution kernels Uh , Uj, and [/§: 

— = U h (Q, Ph, n) x ( 7 ) 

(To dG 6 

Jdsds' J t (s,hj) Uj(s - s' 

J d kUg(k,n,Hs) s d (^Y~ - ^-k,nsj- 


The evolution kernels Uh, Uj and Ug resum the 
large logarithms, ln(C/ 6 ), and explicit expressions can 
be found in Ref. m- The only unknown piece in our 
resummation of logarithms at N 3 LL order is the small 
contribution from the four-loop cusp anomalous dimen¬ 
sion, r™ sf / which we estimate using a Pade approxima¬ 
tion and conservatively vary in our analysis. 

While Eq. (J7| gives the part of the cross section that 
is singular and non-integrable as C —► 0, we also need to 
include the integrable or nonsingular contribution. This 
can be written as 


1 dij ns 
(Tq ~dC~ 


a s (Mns) 
27T 

a s ns) 

27T 


fi(C) 


( 8 ) 


5 (Mns) 


27T 

\fh 

2 

+ 0(a 4 s ). 


/ 2 (C) + /? 0 ln( / ^)/ 1 (C') 
/ 3 (C') + 2/3 0 ln(^)/ 2 (C') 




Q 


fi(C) 


One important property of this shape function is that its 
first moment encodes the leading power correction to our 
cross section. In the MS scheme this moment is given by 

cl kkFc(k). (11) 

Up to the normalization factors shown in Eq. we 
expect approximate universality between the S2i for 
C-parameter and thrust. For the calculation of the cross 
section, the shape function is expanded in a complete ba¬ 
sis of functions obtained by an appropriate infinite-range 
mapping of the Legendre polynomials 133 . with the co¬ 
efficients chosen to maintain the first moment. For fur¬ 
ther details on the implementation of the shape function 
for C-parameter see Ref. [T2]. We remove an O(Aqcd) 
renormalon by using the Rgap scheme uni mu, which in¬ 
troduces a subtraction scale R into our formula, as well 
as the gap parameter A and the perturbative scheme- 
change gap parameter S(R,ns)- Here, S(R,ns) is given 
by a perturbative series in a s (ns) whose mass dimension 
is set by an overall factor of R , and which also contains 
In (ns/R) factors. The convolution with the shape func¬ 
tion now becomes, 

dcr f , o_«(*./*„) 9 d(f / . p\ . 

dc = / d ’’ e 0 "doU-D (12) 

*Fc(p - 37 tA (R,ns)) ■ 

The final component of our cross section is properly 
accounting for hadron-mass effects following Ref. ~HT| . 
Hadron-mass effects induce an additional series of 
large perturbative logarithms which start at NLL, 
c/ ln fc (/rs/AQCD), and also break the exact universality 
between 12/ and S2[ given in Eq. ([ 3 ]). These effects are 
accounted for by including dependence on the transverse 

velocity, r = , /' 1 ,, , in the nonperturbative matrix 

VP±+m 2 H 

elements (here, m# gives the non-zero hadron mass). In 
particular, in the Rgap scheme the first moment of the 
shape function is actually given by 



The functions f\(C), / 2 (C), and / 3 (C) were determined 
in Ref. [12] using the fixed-order results at ©(a 3 ’ 2,3 ) [3] 
HEM EM- The nonsingular cross section d(t ns /d C 
is independent of the renormalization scale p, order by 
order, and therefore we evaluate these pieces at the non¬ 
singular scale /i ns , and vary this scale to estimate higher- 
order perturbative nonsingular corrections. The scale 
variation of p ns will be discussed further in Section [II B 
The full partonic cross section is then given by 


1 do 1 d(t s 1 d(t ns 
apdC = a\)~dC + ~dC 


(9) 


Nonperturbative effects are included by convolving 
Eq. (|9| with a shape function: 


§)*<»■ « 


JdkkF c (k-3TrA(R : ns)) =flf(i?,Ms) (13) 

= 37 t drg c (r)Sli(R,iJ, S ,r). 

Jo 

In the MS scheme the definition accounting for hadron- 
mass effects is the same as Eq. (13), but one sets A = 0 


and removes R as an argument for these parameters. Ac¬ 
counting for both the the MS running due to hadron 
masses and the R-evolution running in the Rgap scheme, 
the evolution of the integrand on the right-hand side of 
Eq. (13) is given by, 


gc(r)Cli(R,n,r) = g c (r ) 


Qis(n) 

a s {fXA) 


1 7i(r) 


l diff 


(RA,R,nA,nu) 


fli(RA,PA,r) 
(14) 
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where i?A and /ta give the initial scales where the func¬ 
tion Qi(Ra, ft a, r) is defined. The perturbative evolution 
kernel A dlff gives the R and g, running for each value of 
r. The function gc encodes the event-shape dependence 
of the hadron-mass effects and 71 gives the solution to 
the one-loop RGE for Hi with hadron masses derived 
in Ref. [33]. Since the two- and three-loop r-dependent 
anomalous dimensions are unknown, we do not treat the 
logs generated by hadron-mass effects to the same level 
of precision. When hadron-nrass effects are accounted 
for we always sum the associated logarithms at NLL. An 
analogous formula to Eq. (14) also holds for the thrust 
parameter fij. 

Combining all of these elements gives the complete 
cross section. Note that we can resum to any order 
up to N 3 LL' and can choose to include or leave out the 
shape function, renormalon subtraction and hadron-mass 
effects. This flexibility allows us to see how the analysis 
changes when we take into account each of these addi¬ 
tional physical considerations and enables us to test how 
robust the fits are to various changes in the theoretical 
treatment. 


B. Profile Functions 


In order to smoothly transition between the nonpertur- 
bative, resummation, and fixed-order regions we make 
use of profile functions g-i(C) for the renormalization 
scales hh, gj{C), gs{C), R(C), and /x ns (C). In the three 
C regions, the requirements on the scales which properly 
deal with large logarithms, nonperturbative effects, and 
the cancellations between singular and nonsingular con¬ 
tributions in the fixed-order region are 


1) nonperturbative: C < 37r Aqcd 

I^h ~ Q, gj ~ V^QCd Q, gs~R~ Aqcd , 

2) resummation: 37 tAqcd CO< 0.75 (15) 

/ C QC 

—, gs~R~ -g- > Aqcd , 

3) fixed-order: C > 0.75 

Mns = Hh = fJ’J = gs = R ~ Q Aqcd • 


In addition we take the fixed-order nonsingular scale 
g „ s ~ gn in the nonperturbative and resummation re¬ 
gions. Our profile functions Hi{C) satisfy these con¬ 
straints, and provide continuous and smooth transitions 
between these C regions. The resummed perturbative 
cross section is independent of 0 ( 1 ) variations in all 
renormalization scales order by order in the logarith¬ 
mic resummation. Therefore the dependence on parame¬ 
ters appearing in the profile functions gets systematically 
smaller as we go to higher orders, and their variation pro¬ 
vides us with a method of assessing perturbative uncer¬ 
tainties. 

For the hard renormalization scale we use the 


C-independent formula 

gn = e-H Q , 


(16) 


where e# is a parameter that we vary from 0.5 to 2 in 
order to account for theory uncertainties. 

The soft scale has different functional dependence in 
the three regions of Eq. (151, and hence depends on the 
following parameters: 


Ms = Hs(C, g 0 ,r s , g H ,t 0 ,t 1 ,t 2 ,t s ). (17) 

Here, go controls the intercept of the soft scale at C = 0, 
t 0 is near the boundary of the purely nonperturbative 
region and t\ controls the end of this transition, where 
the resummation region begins. The transition from non¬ 
perturbative to perturbative is Q dependent, so we use 
the Q-independent parameters no = toQ/(lGeV) and 
n 1 = tiQ/(lGeV). In the resummation region the pa¬ 
rameter r s determines the slope of the soft scale rela¬ 
tive to the canonical resunnnation region scaling, with 
gs = r s gHC/6. The parameter t 2 controls where the 
transition occurs between the resummation and fixed- 
order regions and t s sets the value of C where the renor¬ 
malization scales all become equal as required in the 
fixed-order region. For the jet scale we have the depen¬ 
dence 


MJ = M./(M.ff; Ms(C),ej) , (18) 


where ej is a parameter that is varied in our theory scans 
to slightly modify the natural relation between the scales. 
The exact functional form for gs and [ij in Eqs. © and 
(18) is given in Appendix [A] 

To avoid large logarithms in the soft function subtrac¬ 
tions S, the scale R{C) is chosen to be the same as ns{C) 
in the resummation and fixed-order regions. In the non¬ 
perturbative region we need R{C) < ns(C) to obtain 
an 0(a s ) subtraction that stabilizes the soft function in 
this region (removing unphysical negative dips that ap¬ 
pear in the MS scheme). This introduces an additional 
parameter Rq = R(C = 0). Therefore we have 


R = R(fi S (C),Ro). 


(19) 


The exact functional form for R is also given in Ap¬ 
pendix [A] 

For the nonsingular scale /i ns , we use the variations 

f |[ Vh(C) + hj(C)] , n s = 1 
Mns (C) = < fi H , n s = 0 • ( 20 ) 

[5 [3 (jl h {C) - hj{C)] , n s =- 1 


Here the three choices vary g, ns in a manner that allows 
it to have some independence from gu in the resumma¬ 
tion and nonperturbative regions, while still being equal 
g,H in the fixed-order region (where g,j = g-n)- These 
variations of /i ns probe the higher-order fixed-order un¬ 
certainty in the nonsingular cross section contribution. In 
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parameter 

default value 

range of values 

Mo 

1.1 GeV 

- 

Ro 

0.7 GeV 

- 

n 0 

12 

10 to 16 

n 1 

25 

22 to 28 

t2 

0.67 

0.64 to 0.7 

ts 

0.83 

0.8 to 0.86 

r 3 

2 

1.78 to 2.26 

ej 

0 

-0.5 to 0.5 

ejf 

1 

0.5 to 2.0 

n s 

0 

-1, 0, 1 

j-lCUSp 

1553.06 

- 1553.06 to +4659.18 

S2 

-43.2 

-44.2 to -42.2 

J3 

0 

- 3000 to + 3000 

S3 

0 

- 500 to + 500 

_low 

£2 

0 

-1, 0, 1 

high 

C 2 

0 

-1, 0, 1 

low 

€3 

0 

-1, 0, 1 

high 

e 3 

0 

-1, 0, 1 


TABLE I. C-parameter theory parameters relevant for esti¬ 
mating the theory uncertainty, their default values and range 
of values used for the theory scan during the fit procedure. 
The last four parameters control the statistical errors induced 
by fit functions used in the non-singular terms at O(Qs) (£ 2 ° w 
and £ 2 lgh ) an d 0(a^) (e 3 ° w and eg lsh ) in the region below (e 2 ° w 
and £ 3 ° w ) and above (e !) 1811 and £ 3 Ish ) the shoulder; see Sec. V 
of Ref. pj]. 



FIG. 1. Bands for the profile functions for the renormalization 
scales fiH, hs(C) when varying the profile parameters. 


parameter 

default value 

range of values 

Mo 

1.1 GeV 

- 

Ro 

0.7 GeV 

- 

n 0 

2 

1.5 to 2.5 

ni 

10 

8.5 to 11.5 

t2 

0.25 

0.225 to 0.275 

ts 

0.4 

0.375 to 0.425 

r s 

2 

1.77 to 2.26 

ej 

0 

- 1.5 to 1.5 

eH 

1 

0.5 to 2.0 

n s 

0 

- 1 , 0 , 1 

J3 

0 

- 3000 to + 3000 

S3 

0 

- 500 to + 500 

£ 2 

0 

- 1 , 0 , 1 

£3 

0 

- 1 , 0 , 1 


TABLE II. Thrust theory parameters relevant for estimating 
the theory uncertainty, their default values and range of values 
used for the theory scan during the fit procedure. The last 
two parameters control the statistical errors induced by fit 
functions used in the non-singular terms at 0 ( 0 %) (£ 2 ) and 
0(a 3 ) (£ 3 ); see Sec. E of Ref. [9]. 


vertical arrow on the hard scale indicates the overall 
up/down variation, which causes a variation to all the 
scales. Also shown (as gray dashed lines) are plots of 
the canonical soft scale QC/6 and canonical jet scale 
Q\JC /6. In the resummation region, these correspond 
fairly well with the profile functions, indicating that in 
this region our analysis will avoid large logarithms. As 
discussed in detail in Ref. m, to improve the conver¬ 
gence of the C-parameter cross sections we take r s = 2 as 
our default slope parameter, explaining why our soft and 
jet scales are larger than the canonical values in Fig. [l] 
(by a factor that does not induce further large logs). In 
the analysis of Sec. [Vj we will see how varying each of 
these profile parameters affects our final fit results. 

For the numerical analyses carried out in this work 
we have created within our collaboration two completely 
independent codes. One code within Mathematica [62] 
implements the theoretical expressions exactly as given 
in Ref. |12| , and another code is based on theoretical for¬ 
mulae in Fourier space and realized as a fast Fortran [63] 
code suitable for parallelized runs on computer clusters. 
These two codes agree for the C-parameter distribution 
at the level of 10~ 6 . 


the fixed-order region the variation of /j>h alone precisely 
reproduces the standard fixed-order scale variation. 

The details of the variations of the profile function pa¬ 
rameters used to assess uncertainties are given in Tab. |T| 
The plot in Fig. [l] shows how the scales vary with the 
changes to our C-parameter profile parameters. The 


We will also repeat the thrust fits of Ref. [9], imple¬ 
menting the same type of profile functions used here. 
These profiles have several advantages over those in 
Ref. [9], including a free variable for the slope, a flat 
nonperturbative region, and parameters whose impact 
is much more confined to one of the three regions in 
Eq. (15). For the thrust profiles we redefine r s —> 6r s , 
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which eliminates the factors of 6 in Eqs. © and ( |A4| ). 
This way, the canonical choice of slope is 77 = 1 for both 
C-parameter and thrust. We use r s = 2 as our default for 
thrust as well, again to improve the perturbative conver¬ 
gence, as discussed in Ref. [T5j. The profile parameters 
for thrust and their variations are summarized in Tab. El 
These choices create profiles that are very similar to those 
used in Ref. [TJj. The new fit results for thrust are fully 
compatible with those of Ref. [5] in the resummation re¬ 
gion used for the a s fits. Additionally, they give a better 
description in the nonperturbative region which is out¬ 
side of our fit range. 


C. Hadron-mass effects 


3 7rgc(r)Ui (Ra./A.i) 


Uf (Ra.7A) 



In Ref. mj it was shown that hadron-mass effects, 
apart from breaking the universality properties of the 
leading power correction for various event shapes, also 
induce a nontrivial running. Since these are single loga¬ 
rithms, they start at NLL order. In Ref. [43| the corre¬ 
sponding leading anomalous dimension was determined, 
which yields the NLL resummation of larger logs be¬ 
tween the scales ns and Aqcd for a large set of event 
shapes. The pieces necessary for a higher-order resum¬ 
mation have not yet been computed. One might be wor¬ 
ried that accounting for only the NLL running for Hi in 
an expression as accurate as N 3 LL in cross-section logs 
could be inadequate, or that it could leave significant 
perturbative uncertainties. However, one should recall 
that the hadronic parameter Hi itself is a correction, and 
hence it is valid to account for the related resummation 
with less precision. In this section we show that the Hi 
evolution at NLL order suffices for the precision of our 
N 3 LL' + C>(a 3 ) analysis. Indeed, it turns out that the ef¬ 
fect of the hadron mass running on the fit outcome is very 
small as compared to other uncertainties, and therefore 
can be safely neglected. 

For our C-parameter analysis the implementation of 
hadron mass running effects has been explained at length 
in Ref. T2], and we only summarize here the most rele¬ 
vant aspects needed to understand the fit results. In the 
MS scheme the leading power correction can be written 
as an integral of a universal hadron function, Hi(/r,r), 
common to all event shapes 

^ i ( m ) = Ce f drg e (r)Qi(fj.,r), ( 21 ) 

Jo 

where r is the transverse velocity, e denotes a specific 
event shape, c e is a calculable normalization factor, and 
g e {r) is an event-shape-dependent function encoding the 
hadron-mass effects. The functions g e (r) are known ana¬ 
lytically, and specific examples can be found in Ref. I43| . 
For C-parameter cc = Sir, while for thrust c T = 2. For 
the simple case of the MS scheme the running between 
the initial reference scale no where the universal hadron 
function is specified, and the soft scale ns, is given at 


FIG. 2. Plots of the r dependence of gc(?) 0 .i(Ra, /.ia, r) for 
different values of 9{Ra, Ha)- We normalize to fif (Aa,Pa), 
since it is simply an overall multiplicative factor. 


leading order by 


^1 (Hs,r) = fl l (no,r) 


a s {Vs) 

oisino). 


7i(r) 


( 22 ) 


with 71 (r) = ^ a -ln(l — r 2 ). The corresponding evolution 
formula for the Rgap scheme is considerably more com¬ 
plex, as shown by the form displayed in Eq. (14) above. 


Ensuring that the renormalon is not reintroduced by the 
renormalization group evolution requires an additional 
evolution in the scale R, so A dlff (Ra, R, a*a, r) con¬ 
tains evolution in both the n and R scales. Also here we 
have two reference scales n<\ and Ra to specify the ini¬ 
tial parameter Hi (I?a , nA, v). The full formula for A dlff 
is given in Eq. (67) of Ref. \TZ\ . 

Note that the renormalization group evolution is a 
function of r and takes place independently for each r, 
but the required result for C-parameter or thrust requires 
an integral over r. Due to this integration the functional 
form that we assume for the initial condition Hi ( 7 ', no) 
or Hi(i?Aj g-A, r ) gets entangled with the perturbative 
resummation. 

With current constraints on the r dependence, and 
with the lack of even more precise experimental data to 
probe this issue, a model-independent formulation (like a 
complete functional basis for the r dependence at the ref¬ 
erence scales, Ra and /ta) is not feasible. To implement 
this running we have therefore tested several models for 
the r dependence in Ref. jT2] , and found that generically 
the experimental data is sensitive to the normalization 
which specifies Hi, but not to the detailed form used for 
the r dependence as long as it satisfies several reason¬ 
able constraints. Therefore for the fits performed here 
we simply adopt the default form from Ref. IT2], 

Hi(i? A ,AtA,r) = [a(R A ,nA)fa(r) + b(R A , n^)fb{r)] 2 , 
f a (r) = 3.510e _T3 ^, 


( 23 ) 









f b (r) = 13.585 e ^ - 21.687 e ^. 

This model ensures that P A, J") is always positive 

definite and smoothly goes to zero at the r = 1 endpoint 
(where the ratio of the hadron mass to pt goes to zero). 
The functions / 0i b form an orthonormal basis upon inte¬ 
gration with gc{f)i which yields the following relation: 

FIi(Ra,Pa) = 37t [a(i?A, Pa) 2 + b(R,A, Pa) 2 ], (24) 

which determines the normalization. We also define 

fl(fl A .„ A ) S arct»n ( ^’^ > ) ). (25) 

which was chosen to have an effect orthogonal to the 
more relevant parameter fif (Ra, Pa)- By examining 
our ability to simultaneously measure f2f(i?A> Pa) and 
0(Ra, Pa) we have a means to test for the impact that 
the initial r dependence has on our fits. As we can see 
in Fig. [2j our model captures different behavior for the r 
dependence of fli {Ra , Pa , r) by choosing different values 
of 6(R&,pa). Over the interval r £ [0,1], all the curves 
in this figure are normalized so that they integrate to 1. 


III. EXPERIMENTAL DATA 

Data on the C-parameter cross section are given by 
several experiments for a range of Q from 35 to 207 
GeV. We use data from ALEPH^jwith Q = {91.2, 133, 
161, 172, 183, 189, 200, 206} GeV |M], DELPHI with 
Q = {45, 66, 76, 89.5, 91.2, 93, 133, 161, 172, 183, 
189, 192, 196, 200, 202, 205, 207} GeV [55H58], JADE 
with Q = {35, 44} GeV jMj, L3 with Q = {91.2, 130.1, 
136.1, 161.3, 172.3, 182.8, 188.6, 194.4, 200.2, 206.2} GeV 
[701 HQ, OPAL with Q = {91, 133, 177, 197} GeV [72], 
and SLD with Q = 91.2 GeV [73 1 . As each of these 
datasets is given in binned form, our cross section in 
Eq. @ is integrated over each bin before being com¬ 
pared to the data. The default range on C used in fit¬ 
ting the data is 25 GeV/Q < C < 0.7. A lower limit 
of 25 GeV/Q eliminates the peak region where higher 
nonperturbative moments become important. The 
upper limit is chosen to be 0.7 in order to avoid the far- 
tail region as well as the Sudakov shoulder at C = 0.75. 
Any bin that contains one of the end points of our range 
(C = 25GeV/<2 or 0.7) is included if more than half of 
that bin lies within the range. Using the default range 
and datasets gives a total of 404 bins. As a further test of 
the stability of our analysis, both this C-parameter range 


2 The ALEPH dataset with Q = 91.2 GeV has two systematic 
uncertainties for each bin. The second of these uncertainties 
is treated as correlated while the first one is treated as an un¬ 
correlated uncertainty and simply added in quadrature to the 
statistical uncertainty. 


Q = 35 GeV 




Q = 206 GeV 



FIG. 3. Difference between the default cross section and the 
cross section varying only one parameter. We vary a 3 (mz) 
by ±0.001 (solid red), 2fh by ±0.1 (dashed blue) and 0.2 
by ±0.5 (dotted green). The three plots correspond to three 
different center of mass energies: (a) Q = 35 GeV, (b) Q = 
91.2 GeV, (c) Q = 206 GeV. 


and the selection of datasets is varied in the numerical 
analysis contained in Sec. [V] 

In our fitting procedure, we consider both the sta¬ 
tistical and systematic experimental uncertainties. The 
statistical uncertainties can be treated as independent 
between bins. The systematic experimental uncertain¬ 
ties come from various sources and full documentation of 
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their correlations are not available, so dealing with them 
in our y 2 analysis is more complicated, and we have to use 
a correlation model. For this purpose we follow the LEP 
QCD working group [Ml El and use the minimal overlap 
model. Within one C-parameter dataset, which consists 
of various C-parameter bins at one Q value for one ex¬ 
periment, we take for the bin i, bin j off-diagonal entry 
of the experimental covariance matrix [min(A^ ys , A^ ys )] 2 . 
Here A^ ys are the quoted systematic uncertainties of the 
bins i and j. Within each dataset, this model implies 
a positive correlation of systematic uncertainties. In ad¬ 
dition to this default model choice, we also do the fits 
assuming uncorrelated systematic uncertainties, in order 
to test whether the minimal overlap model introduces 
any bias. See Sec. [VB] for more details on the correlation 
matrix. 


IV. FIT PROCEDURE 

In order to accurately determine both a s (mz) and the 
leading power correction in the same fit, it is important 
to perform a global analysis, that is, simultaneously fit¬ 
ting C-spectra for a wide range of center-of-mass energies 
Q. For each Q , effects on the cross sections induced by 
changes in a s (rriz ) can be partly compensated by changes 
in Q-|, resulting in a fairly strong degeneracy. This is re¬ 
solved by the global fit, just as in the thrust analysis of 
Ref. [9]. Fig. [3] shows the difference between the theoret¬ 
ical prediction for the cross section at three different val¬ 
ues of Q , when a s (mz) or Hi are varied by ±0.001 and 
±0.05 GeV, respectively. It is clear that the potential 
degeneracy in these parameters is broken by having data 
at multiple Q values. In Fig. [3] we also vary the higher- 
order power correction parameter fi^ , which clearly has a 
much smaller effect than the dominant power correction 
parameter Hi. 

To carry out a fit to the experimental data we fix 
the profile and theory parameters to the values shown 
in Tab. [Tj The default values are used for our primary 
theory cross section. We integrate the resulting theo¬ 
retical distribution over the same C-parameter bins as 
those available experimentally, and construct a y 2 func¬ 
tion with the uncorrelated statistical experimental uncer¬ 
tainties and correlated systematic uncertainties. This y 2 
is a function of a s (mz) and Hi, and is very accurately 
described by a quadratic near its global minimum, which 
therefore determines the central values and experimental 
uncertainties. The value of fR and its associated uncer¬ 
tainties encode the dominant hadronization effect as well 
as the dominant residual uncertainty from hadronization. 

To obtain the perturbative theoretical uncertainty we 
consider the range of values shown for the theory parame¬ 
ters in Tab.[j] Treating each of these as a flat distribution, 
we randomly generate values for each of these parameters 
and then repeat the fit described above with the new y 2 
function. This random sampling and fit is then repeated 
500 times. We then construct the minimum ellipse that 


fully contains all 500 of the central-fit values by first cre¬ 
ating the convex envelope that contains all of these points 
within it. Then, we find the equation for the ellipse that 
best fits the points on the envelope, with the additional 
restrictions that all values lie within the ellipse and its 
center is the average of the maximum and minimum val¬ 
ues in each direction. This ellipse determines the per¬ 
turbative theoretical uncertainty, which turns out to be 
the dominant uncertainty in our fit results. In our final 
results the perturbative and experimental uncertainties 
are added in quadrature. This procedure is similar to 
that discussed in the Appendix of Ref. uni. 


V. RESULTS 


In this section we discuss the results from our global 
analysis. We split the presentation into several subsec¬ 
tions. In Sec. |V A] we discuss the impact that resumma¬ 
tion and the inclusion of power corrections have on the fit 
results. In Sec. |V B| we present the analysis which yields 
the perturbative uncertainty in detail, cross-checking our 
method by analyzing the order-by-order convergence. We 
also analyze the impact of removing the renormalon. In 
Sec. |V C| we discuss the experimental uncertainties ob¬ 
tained from the fit. Section VD| discusses the impact 
that varying the theory parameters one by one has on 
the best-fit points, allowing us to determine which pa¬ 
rameters dominate the theoretical uncertainty. The im¬ 
pact of hadron-mass resummation is discussed in detail 
in Sec. |VE| We examine the effects of changing the de¬ 
fault dataset in Sec. IV FI The final fit results are collected 
in Sec. V G When indicating the perturbative precision, 
and whether or not the power correction Hi is included 
and at what level of precision, we use the following nota¬ 
tion: 


0(a k s ) 

N fc LL '+0{a k ) 
N fe LL'±0(a*)±fii 
N fc LL'±C?(aJ)+fii(.R, n) 
~N k LL'+0{a k s )+n 1 (R, n,r) 


fixed order up to 0(a k ) 
perturbative resummation 
MS scheme for Hi 
Rgap scheme for Hi 
Rgap scheme with 
hadron masses for Hi. 


A. Impact of Resummation and Power Corrections 

In Fig. [4] we show a s (mz) extracted from fits to the 
tail of the C-parameter distribution including sequential 
improvements to the treatment of perturbative and non- 
perturbative components of our code, using the highest 
perturbative accuracy at each stage. The sequence from 
left to right shows the fit results using: 0(a 3 s ) fixed-order 
results only, adding N 3 LL resummation, adding the H x 
power correction, adding renormalon subtractions and 
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a s {mz ) from global C-parameter tail fits 




r 


_ 

0.135 


1 

Below error bars & ± -» perturbative error 

- 

a s (m z ) 


0 ( 0 ^) fixed-order 

All errors: ff s ( m z) = 0.1123 ±0.0015 

- 



1 0.1317 ±0.0052 



0.130 


1 


— 

0.125 



i 

i 


- 



1 + N 3 LL’ summation + Power Collection 

- 

0.120 


0.1219 ±0.0028 

0.1117 ±0.0016 

— 



1 

+ R-scheme + hadron mass effects - 



1 

0.1123 ±0.0014 0 . 1119 ± 0.0013 

- 

0.115 


1 


- 



1 

| } i 

- 

0.110 


_L 


- 


FIG. 4. The evolution of the value of a a (mz) adding components of the calculation. An additional ~ 8% uncertainty from not 
including power corrections is not included in the two left points. 


using the Rgap power correction parameter Qi(Ra, Ha), 
and adding hadron-mass effects. These same results to¬ 
gether with the corresponding x 2 /dof are also collected 
The fit with only fixed-order 0(al) results 


III 


in Tab. 

has a relatively large y 2 /dof and also its central value 
has the largest value of a s (mz)- Including the resum¬ 
mation of large logarithms decreases the central a s (mz) 
by 8% and also decreases the perturbative uncertainty 
by ~ 50%. Due to this smaller perturbative uncertainty 
it becomes clear that the theoretical cross section has 
a different slope than the data, which can be seen, for 
example, at Q = mz for 0.27 < C < 0.35. This leads 
to the increase in the y 2 /dof for the “N 3 LL' no power 
corr.” fit, and makes it quite obvious that power correc¬ 
tions are needed. When the power correction parameter 
fli is included in the fit, shown by the third entry in 
Tab. |III| and the result just to the right of the vertical 
dashed line in Fig. [4j the y 2 /dof becomes 1.004 and this 
issue is resolved. Furthermore, a reduction by ~ 50% 
is achieved for the perturbative uncertainty in a s (rnz ). 
This reduction makes sense since some of the perturba¬ 
tive uncertainty of the cross section is now absorbed in 
Q-i, and a much better fit is achieved for any of the vari¬ 
ations associated to estimating higher-order corrections. 
The addition of also caused the fit value of a s (mz) to 
drop by another 8%, consistent with our expectations for 
the impact of power corrections and the estimate made in 
Ref. [T2]. Note that the error bars of the first two purely 
perturbative determinations, shown at the left-hand side 
of the vertical thick dashed line in Fig. [I] and in the last 
two entries in Tab. Ill do not include the ~ 8% uncer¬ 


tainties associated with the lack of power corrections. 

The remaining corrections we consider are the use of 
the R-scheme for fli which includes the renormalon sub¬ 


tractions, and the inclusion of the log-resummation ef¬ 
fects associated to the hadron-mass effects. Both of these 
corrections have a fairly small impact on the determi¬ 
nation of a s (mz), shifting the central value by +0.5% 
and — 0.3% respectively. Since adding the — 0.3% shift 
from the hadron mass corrections in quadrature with the 
~ 1.2% perturbative uncertainty does not change the 
overall uncertainty we will use the R-scheme determi¬ 
nation for our main result. This avoids the need to fully 
discuss the extra fit parameter 9(Ra,Ha) that appears 
when hadron masses are included. Further discussion 
of the experimental uncertainties and the perturbative 
uncertainty from the random scan are given below in 
Secs. VB and VD and a more detailed discussion of 
the impact of hadron-mass resummation is given below 
in Sec. IV El 

The values of fli obtained from the fits discussed above 
can be directly compared to the fli power correction ob¬ 
tained from the thrust distribution. Values for flj from 
the C-parameter fits are given below in Secs. |V Bl and |V D| 
and the comparison with thrust is considered in Sec.|VII| 


B. Perturbative Uncertainty from the Scan 

To examine the robustness of our method of determin¬ 
ing the perturbative uncertainty by the random scan, we 
consider the convergence and overlap of the results at dif¬ 
ferent perturbative orders. Figure [5lshows the spread of 
best-fit values at NLL', N 2 LL' and N^LL'. The upper left 
panel, Fig. |5(a)[ shows results from fits performed in the 
Rgap scheme, which implements a renormalon subtrac¬ 
tion for f2 1; and the upper right-panel, Fig. |5(b)| shows 
results in the MS scheme without renormalon subtrac- 
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a s( m z) 

(a) 



a s {m z ) 

(b) 


p—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—i—r- 


1.4- 


1.3- 


r_ 

dof 


1 . 2 - 


1.1 - 


N 3 LL' + 0(a 3 ) + n,(R,p) 

n 2 ll' + o(ff 2 ) + ni(R,^) 

NLL' + Ola's) + fli(R,p) 

Full Results 

T\ 

\ 


• \ 




■ * i. ‘ ‘ T 5 - ■ 

V.V 



i .... i 


i .... i 


i .... i 


0.1000 0.1025 0.1050 0.1075 0.1100 0.1125 0.1150 

a s (m z ) 

(c) 



a s (m z ) 

(d) 


FIG. 5. The first two panels show the distribution of best-fit points in the a s (rnz)- 20i and a 3 (jnz)~ 2f2i planes. Panel (a) 
shows results including perturbation theory, resummation of large logs, the soft nonperturbative function and ffi defined in 
the Rgap scheme with renormalon subtractions. Panel (b) shows the results as in panel (a), but with fR defined in the MS 
scheme without renormalon subtractions. In both panels the dashed lines corresponds to an ellipse fit to the contour of the 
best-fit points to determine the theoretical uncertainty. The respective total (experimental + theoretical) 39% CL standard 
uncertainty ellipses are displayed (solid lines), which correspond to 1-cr (68% CL) for either one-dimensional projection. The 
big points represent the central values in the random scan for a s (mz) and 2fii. Likewise, the two panels at the bottom 
show the distribution of best-fit points in the a 3 {m z )-X 2 /dof plane. Panel (c) shows the x 2 /dof values of the points given in 
panel (a), whereas panel (b) shows the x 2 /dof values of the points given in panel (b). 


tions. Each point in the plot represents the outcome of orders in perturbation theory. Not unexpectedly, fits in 
a single fit, and different colors correspond to different the Rgap scheme show generally smaller theory uncer- 












12 




a s (m z ) 

X 2 /dof 

n 3 ll' 

+ hadron 

0.1119(13)(06) 

0.991 

n 3 ll' 

with fli (R, fi) 

0.1123(14)(06) 

0.988 

n 3 ll' 

with Sli 

0.1117(16)(06) 

1.004 

n 3 ll' 

no power corr. 

0.1219(28)(02) 

2.091 

O(cts) fixed order 
no power corr. 

0.1317(52)(03) 

1.486 


TABLE III. Comparison of C-parameter tail fit results for 
analyses when we add various components of the theoreti¬ 
cal result (from the bottom to top). The first parentheses 
gives the theory uncertainty, and the second is the experi- 


mental and hadronic uncertainties added in quadrature for 
the first three rows, and experimental uncertainty for the last 
two rows. 

order 

a s (mz) (with Hi) 

a 3 {mz) (with Qi(Ra, ^a)) 

NLL' 

0.1071(60)(05) 

0.1059(62) (05) 

n 2 ll' 

0.1102(32)(06) 

0.1100(33)(06) 

N 3 LL' (full) 

0.1117(16)(06) 

0.1123(14)(06) 


TABLE IV. Central values for a s (mz) at various orders with 
theory uncertainties from the parameter scan (first value 
in parentheses), and experimental and hadronic uncertainty 
added in quadrature (second value in parentheses). The bold 
N 3 LL' value is our final result. 


tainties. 

In order to estimate correlations induced by theoreti¬ 
cal uncertainties, each ellipse in the a s -2I7i plane is con¬ 


structed following the procedure discussed in Sec. IV 


Each theory ellipse constructed in this manner is inter¬ 
preted as an estimate for the i-a theoretical uncertainty 
ellipse for each individual parameter (39% confidence for 
the two parameters), and is represented by a dashed el¬ 
lipse in Fig. [5] The solid lines represent the combined 
(theoretical plus experimental) standard uncertainty el¬ 
lipses at 39% confidence for two parameters, obtained by 
adding the theoretical and experimental error matrices 
from the individual ellipses, where the experimental el¬ 
lipse corresponds to Ay 2 = 1. Figure [ 5 ] clearly shows 
a substantial reduction of the perturbative uncertainties 
when increasing the resummation accuracy, and given 
that they are 39% confidence regions for two parameters, 
also show good overlap between the results at different 
orders. 

The results for a s (mz) and 17 1 from the theory scan at 
different perturbative orders are collected in Tabs |IV| and 
|Vj Central values here are determined from the average 
of the maximal and minimal values of the theory scan, 
and are very close to the central values obtained when 
running with our default parameters. The quoted per¬ 
turbative uncertainties are one-parameter uncertainties. 

In Tab. Ill above we also present a s (mz ) results with 
no power corrections and either using resummation or 
fixed-order perturbative results. Without power correc- 


order fti [GeV] 

NLL' 0.533(154)(18) 

N 2 LL' 0.443(119) (19) 

N 3 LL' (full) 0.384(91)(20) 


»i(Ra,Ma) [GeV] 
0.582(134)(16) 
0.457(83)(19) 

0.421(60)(20) 


TABLE V. Central values for at the reference scales 
Ra = A* a = 2 GeV and for f2i and at various orders. The 
parentheses show first the theory uncertainties from the pa¬ 
rameter scan, and second the experimental plus the uncer¬ 
tainty due to the imprecise determination of a s (added in 
quadrature). The bold N 3 LL' value is our final result. 


tions there is no fit for Hi, so we take the central value 
to be the average of the maximum and minimum value 
of a s (niz) that comes from our parameter scan. Our 
estimate of the uncertainty is given by the difference be¬ 
tween our result and the maximum fit value. For the 
fixed-order case, since there is only one renormalization 
scale, we know that the uncertainties from our parame¬ 
ter variation for e#, s^, eif™ and £ 3 ™ are uncorrelated. 
So, we take the fit value for ot s (rriz) with the default 
parameters as our result and add the uncertainties from 
variations of these parameter in quadrature to give the 
total uncertainty. 


An additional attractive result of our fits is that the ex¬ 
perimental data is better described when increasing the 
order of the resummation and fixed-order terms. This 
can be seen by looking at the minimal % 2 /dof values 
for t he be st-fi t poi nts, which are shown in Fig. [5] In 
Figs. 5(c) and 5(d) we show the distribution of x% n /dof 


values for the various a s (mz) best-fit points. Figure |5(c) 
displays the results in the Rgap scheme, whereas Fig. 15 (d 


shows the results in the MS scheme. In both cases we 
find that the Xmin values systematically decrease with 
increasing perturbative order. The highest-order analy¬ 
sis in the MS scheme leads to Xmin/dof values around 
unity and thus provides an adequate description of the 
whole dataset, however one also observes that account¬ 
ing for the renormalon subtraction in the Rgap scheme 
leads to a substantially improved theoretical description 
having Xmin/d°f values below unity essentially for all 
points in the random scan. Computing the average of 
the Xmin values we find at N 3 LL' order for the Rgap 
and MS schemes 0.988 and 1.004, respectively (where the 
spread of values is smaller in the Rgap scheme). Likewise 
for N 2 LL' we find 1.00 and 1.02, and for NLL' we find 
1.09 and 1.14. These results show the excellent descrip¬ 
tion of the experimental data for various center-of-mass 
energies. They also validate the smaller theoretical un¬ 
certainties obtained for a s and fR at N 2 LL' and N 3 LL' 
orders in the Rgap scheme. 
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o's(mz) 


T/tot 

V ij 


/ 2.35(5) -lO" 6 -0.51(5) • 10" 4 GeV \ 

\ - 0.51(5) • 10 -4 GeV 1.61(5) ■ 10" 2 GeV 2 J ' 


Note that the theoretical uncertainties dominate by a sig¬ 
nificant amount. The experimental correlation coefficient 
is significant and reads 

Pan = -0.93(15). (27) 


The theory correlation coefficient is small, P^q° = 
— 0.16, and since these uncertainties dominate it reduces 
the correlation coefficient for the total uncertainty to 


PaO = _ 0 . 26 ( 2 ). 


(28) 


In both Eqs. (27) and (28) the numbers in parentheses 


indicate a ± range that captures all values obtained from 
the theory scan. The correlation exhibited by the green 
dotted experimental uncertainty ellipse in Fig. [6] is given 
by the line describing the semimajor axis 


FIG. 6. Experimental Ay 2 = 1 standard uncertainty ellipse 
(dotted green) at N 3 LL/ accuracy with renormalon subtrac¬ 
tions, in the a s - 2fh plane. The dashed blue ellipse represents 
the theory uncertainty which is obtained by fitting an ellipse 
to the contour of the distribution of the best-fit points. This 
ellipse should be interpreted as the 1-a theory uncertainty 
for one parameter (39% confidence for two parameters). The 
solid red ellipse represents the total (combined experimental 
and perturbative) uncertainty ellipse. 


C. Experimental Fit Uncertainty 

Next we discuss in more detail the experimental un¬ 
certainty in a s (mz) and the hadronization parameter fix 
as well as the combination with the perturbative uncer¬ 
tainty done to obtain the total uncertainty. 

Results are depicted in Fig. [6] for our highest order 
fit including resummation, power corrections and renor¬ 
malon subtractions. The inner green dotted ellipse, 
blue dashed ellipse, and solid red ellipse represent the 
Ay 2 = 1 uncertainty ellipses for the experimental, the¬ 
oretical, and combined theoretical and experimental un¬ 
certainties respectively. These ellipses correspond to the 
one-dimensional projection of the uncertainties onto ei¬ 
ther a s (mz ) or fix (39% confidence ellipse for two param¬ 
eters). The correlation matrix of the experimental, the¬ 
ory, and total uncertainty ellipses are (for i,j = a s , 2 fix), 


fix 

30.84 GeV 


= 0.1257 -a s {m z ). 


(29) 


Note that extrapolating this correlation to the extreme 
case where we neglect the nonperturbative corrections 
(fix = 0) gives a s (mz) —> 0.1257 which is consistent 
with the 0.1219 ± 0.0028 result of our fit without power 
corrections in Tab. m 


From V i ® xp in Eq. (|26|) it is possible to extract the 


experimental uncertainty for a s (mz) and the uncertainty 
due to the imprecise determination of fix, 


tin = 0.0002, 

= cr as || = 0.0006, (30) 

and to extract the experimental uncertainty for fix and 
its uncertainty due to the imprecise determination of 
ot s {m z ), 


°nT = o-Ox \A - Pin = °- 014 GeV - 

= <r ai \p a n\ = 0.037 GeV. (31) 

The projections of the outer solid ellipse in Fig. [6] show 
the total uncertainty in our final one-parameter results 
obtained from ld* ot , which are quoted below in Eq. (f34j). 


D. Individual Theory Scan Errors 


Vij — 


G a s 2CTa a Cr ni Pan 

2o’a a cr OiPan 4 (Tq ! 


(26) 


l^ xp = ( 

' 4.18(52) • 10” 7 

-0.24(5) • 10 -4 


v - 0.24(5) • 10~ 4 GeV 1.60(47) • 10" 4 ( 

VV eo = ( 

f 1.93 • 10~ 6 

— 0.27 ■ 10 -4 GeV 

v \ 

, - 0.27 • 10~ 4 GeV 

1.45 • 10” 2 GeV 2 


To gain further insight into our theoretical precision 
and in order to estimate the dominant source for theory 
uncertainty from missing higher-order terms, we look at 
the size of the theory uncertainties caused by the indi¬ 
vidual variation of each one of the theory parameters 
included in our random scan. In Fig. [7] two bar charts 
are shown with these results for a s (mz) (left panel) and 
£2i(i?A,P a) (right panel) for fits corresponding to our 
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-0.0002 0 0.0002 0.0004 0.0006 


Aa's(mz) 



-0.03 -0.02 -0.01 0 0.01 

MMGeV] 


FIG. 7. Variations of the best-fit values for a s (mz) and fli from up (dark blue) and down (light green) variations for the 
theory parameters according to Tab. [I] We do not display those parameters which do not affect the fit region (e!) 18 * 1 , £ 3 lgl \ M 0 j 

Ro, no). 


best theoretical setup (with N 3 LI/ accuracy and in the 
Rgap scheme). The dark blue bars correspond to the 
result of the fit with an upward variation of the given 
parameter from Tab. [TJ while the light green bars corre¬ 
spond to the fit result from the downward variation in 
Tab. [T] Here we vary a single parameter keeping the rest 
fixed at their default values. We do not show parameters 
that have a negligibly small impact in the fit region, e.g. 
e* 2 lsh and eg 1 ® 11 , which only have an effect on the cross 
section to the right of the shoulder, or no , which only 
affects the cross section in the nonperturbative region. 

We see that the dominant theory uncertainties are re¬ 
lated to variations of the profile functions (ejj, r s , ej, £ 2 ), 
where e# is the largest source of uncertainty, and is par¬ 
ticularly dominant for fix. The second most important 
uncertainty comes from r s for a s and G for Q-[. and ej 
also has a significant effect on both parameters. 

As expected, the parameters associated to the tran¬ 
sitions on the sides of our fit region, n\ and t s , hardly 
matter. The renormalization scale parameter n s for the 
nonsingular partonic distribution d<j ns /dC also causes a 
very small uncertainty since the nonsingular terms are 
always dominated by the singular terms in our fit region. 
The uncertainties related to the numerical uncertainties 
of the perturbative constants (s^, s 3 , jo) as well as the 
numerical uncertainties in the extraction of the nonsin¬ 
gular distribution for small C values, (e 2 ° w > e 3 ° w ) are - 
with the possible exception of jo - much smaller and do 
not play an important role. The uncertainty related to 


the unknown 4-loop contribution to the cusp anomalous 
dimension is always negligible. Adding quadratically the 
symmetrized individual uncertainties shown in Fig. [7j we 
find 0.0007 for a s and 0.05 GeV for Oi. This is about 
one half of the theoretical uncertainty we have obtained 
by the theory parameter scan for a s (or five sixths for 
f2i), demonstrating that incorporating correlated vari¬ 
ations through the theory parameter scan represents a 
more realistic method to estimate the theory uncertainty. 


E. Effects of $7i hadron-mass resummation 


The fit results presented in the previous two sections 
ignored the small hadron-mass effects. These effects are 
analyzed in greater detail in this section. We again per¬ 
form 500 fits for a theory setup which includes N 3 LL' ac¬ 
curacy and a power correction in the Rgap scheme, but 
this time it also includes hadron-mass-induced running. 

Since the impact of hadron-mass effects is small, one 
finds that the experimental data in the tail of the dis¬ 
tribution is not accurate enough to fit for 0 (Ra,Ma) in 


Eq. (25), in addition to a s (niz) and Qi{Ra, Ha)- This is 
especially true because it enters as a small modification to 
the power correction, which by itself is not the dominant 
term. Indeed, fitting for a(f?A,MA) and 5(1?a,Ma) as de¬ 
fined in Eq. ( |23| ) gives a strongly correlated determina¬ 
tion of these two parameters. The dominant hadronic pa¬ 
rameter flf (i?A 5 A*a)j which governs the normalization, 
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ff s ( m z ) 

FIG. 8. Comparison of fits to the C-parameter tail distribu¬ 
tion with theory prediction which include/ignore hadron-mass 
effects (in red/blue). Although a direct comparison of a 3 val¬ 
ues is possible, one has to keep in mind that fli(/XA, Ra) has 
a different meaning once hadron mass running effects are in¬ 
cluded. 


is still as accurately determined from data as the fli in 
Tab.|V| However, the orthogonal parameter 9(Ra,^a) is 
only determined with very large statistical uncertainties. 
As discussed in Ref. [T2j, the specific value of 9(Ra,^a) 
has a very small impact on the cross section, which is 
consistent with the inability to accurately fit for it. 

The results of our fit including hadron-mass effects are 


a s (m z ) = 0.1119 ± 0.0006 exp+had ± 0.0013 pert , (32) 

Ma) = 0.411 ± 0.018 exp -t-a s ± 0.052 per t GeV. 


Note that the meaning of H 1 (i?A,M a) here is different 
from the case in which hadron-mass running effects are 
ignored because there are extra evolution effects needed 
to translate this value to that used in the cross section 
at a given value of C, compared to the no-hadron-mass 
case. 

In Fig. [8] we compare the outcome of the 500 fits at 
N 3 LL' in the Rgap scheme. Results with hadron-mass ef¬ 
fects give the red ellipse on the left, and without hadron- 
mass effects give the blue ellipse on the right. (The latter 
ellipse is the same as the one discussed above in Sec. V B ) 
The effects of hadron masses on a s {mz) are to decrease 
its central value by 0.3% and reduce the percent pertur¬ 
bative uncertainty by 0.1%. Given that the total pertur¬ 
bative uncertainties are 1.2%, these effects are not sta¬ 
tistically significant. When studying the effect on Hi 
one has to keep in mind that its meaning changes when 
hadron-mass effects are included. Ignoring this fact we 
observe that hadron masses shift the central value down¬ 


wards by 2.4%, and reduce the percent theoretical un¬ 
certainty by 1.6%. Again, given that the perturbative 
uncertainty for Hi is 14%, this shift is not significant. 


Since the theory uncertainties become slightly smaller 
when hadron-mass effects are incorporated, one could use 
this setup as our default. However we take a more con¬ 
servative approach and consider the 0.3% shift on the 
central value as an additional source of uncertainty, to 
be added in quadrature to the hadronization uncertainty 
already discussed in Sec. VB This increases the value 
of the hadronization uncertainty from 0.0006 to 0.0007, 
and does not affect the total a s uncertainty. The main 
reason we adopt this more conservative approach is that, 
while well motivated, the ansatz that we take in Eq. (23) 
is not model independent. We believe that this ansatz 
serves as a good estimate of what the numerical effect of 
hadron masses are, but should likely not be used for the 
central fit until further theoretical insight on the form of 
Hi(r) is gained. We do not add an additional uncertainty 
to Hi since hadron-mass effects change its meaning and 
uncertainties for Hi are large enough that these effects 
are negligible. 


In App. [B] we also consider fits performed using the 
Rgap scheme with C-parameter gap subtractions, rather 
than our default Rgap scheme with thrust gap subtrac¬ 
tions. The two results are fully compatible. As discussed 
in Ref. [12] the thrust gap subtractions give better per¬ 
turbative convergence, and hence are used for our default 
cross section. 


F. Dataset dependence 


In this section we discuss how much our results de¬ 
pend on the dataset choice. Our default global dataset 
accounts for all experimental bins for Q > 35 GeV in the 
intervals [C ^in,C max [ = [25/Q,0.7], (more details are 
given in Sec. uni) The upper limit in this range is moti¬ 
vated by the fact that we do not want to include data too 
close to the shoulder, since we do not anticipate having 
the optimal theoretical description of this region. The 
lower limit avoids including data too close to the non- 
perturbative region, which is near the cross section peak 
for Q = mz , since we by default only include the leading 
power correction Hj in the OPE of the shape function. 
To consider the impact of this dataset choice we can vary 
the upper and lower limits used to select the data. 

In Fig. [9] the best fits and the respective total exper¬ 
imental + theory 68% CL uncertainty ellipses (for two 
parameters) are shown for global datasets based on dif¬ 
ferent choices of data ranges. The result for our default 
global dataset is given in red, with a thicker, dashed el¬ 
lipse. In the caption of Fig. [9] the data ranges and the 
number of bins are specified for each one of the plotted 
ellipses. 

Interestingly all uncertainty ellipses have very similar 
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a s (m z ) 

FIG. 9. Global fit results for different choices of 
dataset, using our best theory setup at N 3 LL / with 
power corrections in the Rgap scheme. Consider¬ 
ing the central values from left to right, the datasets 
read [ Cmin, Umax ] # of bins : [29/Q, 0 . 7 ) 371 , [ 22 /Q, 0 . 75 ] 453 , 
[23/Q, 0.7]4i7, [0.24,0.75)403, [24/Q, 0.7] 40 9, [25/Q, 0.7 ] 404 
(default), [25/Q, 0 . 6 ) 322 , [ 25/Q, 0.75] 43 o, [27/Q, 0.7] 38 6, 
[25/Q, 0 . 65 ) 349 , [22/Q, 0 . 7 ) 427 . We accept bins which are 
at least 50% inside these fit regions. The ellipses correspond 
to total 1-cr uncertainties (experimental + theory) for two 
variables (a s and Hi), which are suitable for a direct compar¬ 
ison of the outcome of two-parameter fits. The center of the 
ellipses are also shown. 


correlation and are lined up approximately along the line 


f 'll 

41.26 GeV 


0.1221 — a s (mz) ■ 


(33) 


As expected, the results of our fits depend only weakly on 
the C range and the size of the global datasets, as shown 
in Fig. [9] The size and tilt of the total uncertainty el¬ 
lipses is very similar for all datasets (with the exception of 
[22/Q,0.7], which clearly includes too much peak data). 
Since the centers and the sizes of the uncertainty ellipses 
are fully statistically compatible at the 1-cr level, this 
indicates that our theory uncertainty estimate at N 3 LL' 
really reflects the accuracy at which we are capable of de¬ 
scribing the different regions of the spectrum. Therefore 
a possible additional uncertainty that one could consider 
due to the arbitrariness of the dataset choice is actually 
already represented in our final uncertainty estimates. 


G. Final Results 

As our final result for a s (mz ) and fli, obtained at 
N 3 LL' order in the Rgap scheme for Pa), we get 

(34) 


Cd a 
o-dC 



FIG. 10. C-parameter distribution at N 3 LL' order for Q = 
rriz showing the fit result for the values for a s (mz) and fii. 
The blue band corresponds to the theory uncertainty as de¬ 
scribed in Sec. |V B| Experimental data is also shown. 


± 0.0007hadr ± 0.0014 pert , 
^i(Ra,Ma) = 0.421 ± 0.007exp 

± 0.019 aa ( mz ) ± 0.060 P ert GeV, 

where Ra = Pa = 2 GeV and we quote individual 1-cr 
uncertainties for each parameter. Here x 2 /doi = 0.99. 
Equation (34) is the main result of this work. 


Equation (34) accounts for the effect of hadron mass 


running through an additional (essentially negligible) un¬ 
certainty. Also, it neglects QED and finite bottom-mass 
corrections, which were found to be small effects in the 
corresponding thrust analysis in Ref. [9] . 

Given that we treat the correlation of the system¬ 
atic experimental uncertainties in the minimal overlap 
model, it is useful to examine the results obtained when 
assuming that all systematic experimental uncertain¬ 
ties are uncorrelated. At N 3 LL' order in the Rgap 
scheme the results that are analogous to Eq. (34) read 


a s (m z ) = 0.1123±0.0002 exp ±0.0007 had r±0.0012 pert and 
Hi(R a , Pa) = 0.412 ± 0.007 exp ±0.022 as ±0.061 per t GeV 

with a combined correlation coefficient of p/Qp al = 
— 0.091. The results are compatible with Eq. (34), in¬ 


dicating that the ignorance of the precise correlation of 
the systematic experimental uncertainties barely affects 
the outcome of the fit. 

In Fig. [lO] we show the theoretical fit for the 
C-parameter distribution in the tail region, at a center- 
of-mass energy corresponding to the Z- pole. We use the 
best-fit values given in Eq. (34). The band corresponds to 


the perturbative uncertainty as determined by the scan. 
The fit result is shown in comparison with experimental 
data from DELPHI, ALEPH, OPAL, L3 and SLD. Good 
agreement is observed for this spectrum, as well as for 
spectra at other center of mass values. 


a s (m z ) = 0.1123 ± 0.0002 exp 
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0-dC 



FIG. 11. C-parameter distribution below the fit region, shown 
at N 3 LL/ order for Q = mz using the best-fit values for 
a s (rnz) and fii. Again the blue band corresponds to the 
theory uncertainty and error bars are used for experimental 
data. 


1 dcr 
crdC 



FIG. 12. C-parameter distribution above the fit range, shown 
at N 3 LL/ order for Q = mz using the best-fit values for 
a s (mz) and Hi. Again the blue band corresponds to the the¬ 
ory uncertainty and the error bars are used for experimental 
data. 


VI. PEAK AND FAR TAIL PREDICTIONS 


in the peak region are shown in Fig. |11| The nice agree¬ 
ment within theoretical uncertainties (blue band) with 
the precise data from DELPHI, ALEPH, OPAL, L3, and 
SLD indicates that the value of Hi obtained from the fit 
to the tail region is the dominant nonperturbative effect 
in the peak. The small deviations between the theory 
band and the experimental data can be explained due to 
the fact that the peak is also sensitive to higher-order 
power corrections D(? >2 , which have not been tuned to 
reproduc e th e peak data in our analysis. 

In Fig. 12 we compare predictions from our full N 3 LL' 
code in the Rgap scheme to the accurate DELPHI, 
ALEPH, L3, and SLD data at Q = mz in the far-tail re¬ 
gion]^] We find excellent agreement with the data within 
the theoretical uncertainties (blue band). The key fea¬ 
ture of our theoretical prediction that matters most in 
the far-tail region is the merging of the renormalization 
scales toward /i.g = /ij = p h at C ~ 0.75 in the profile 
functions. This is a necessary condition for the cancel¬ 
lations between singular and nonsingular terms in the 
cross section to occur above the shoulder region]^] At 
Q = mz the theoretical cross section presented here ob¬ 
tains accurate predictions in the region both below and 
above the shoulder that agree with the data. Our analy¬ 
sis does not include the full 0 (a k Aqcd/Q) power correc¬ 
tions (for k < 4), since they are not part of our master 
formula. Nevertheless, and in analogy with what was 
found in the case of thrust, agreement with the experi¬ 
mental data seems to indicate that these missing power 
corrections may be smaller than naively expected. 


VII. UNIVERSALITY AND COMPARISON TO 
THRUST 

An additional prediction of our theoretical formal¬ 
ism is the universality of Hi between the thrust and 
C-parameter event shapes. Therefore, a nontrivial test 
of our formalism can be made by comparing our result 
for Hi with the determination from the earlier fits of the 
thrust tail distributions in Ref. |9j and the first moment 
of the thrust distribution in Ref. EDI- 

Since we have updated our profiles for thrust, it is ex¬ 
pected that the outcome of the a s and Hi determination 
is slightly (within theoretical uncertainties) different from 
that of Ref. [3] ■ We also have updated our code to match 


Even though our fits were performed in the resumma¬ 
tion region which is dominated by tail data, our theoret¬ 
ical results also apply for the peak and far-tail regions. 
As an additional validation for the results of our global 
analysis in the tail region, we use the best-fit values ob¬ 
tained for a s and fli to make predictions in the peak and 
the far-tail regions where the corresponding data was not 
included in the fit. 

Predictions from our full N 3 LL' code in the Rgap 
scheme for the C-parameter cross section at the Z -pole 


! The OPAL data was excluded from the plot because its bins are 
rather coarse in this region, making it a bad approximation of 
the differential cross section. 

4 It is worth mentioning that in the far-tail region we employ the 
MS scheme for f2i, since the subtractions implemented in the 
Rgap scheme clash with the partonic shoulder singularity, re¬ 
sulting in an unnatural behavior of the cross section around 
C = 0.75. The transition between the Rgap and MS schemes 
is performed smoothly, by means of a hybrid scheme which in¬ 
terpolates between the two in a continuous way. This hybrid 
scheme has been discussed at length in Ref. m- 
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FIG. 13. Comparison of determinations of a s (mz) and fir 
with the corresponding total 1 -a uncertainty ellipses. As an 
illustration we display the determination of fif obtained from 
fits to the C-parameter distribution (green), which is clearly 
different from obtained from thrust fits (blue), and the de¬ 
termination of fli as obtained from C-parameter distribution 
fits (red). All fits have been performed with N 3 LL' theoretical 
predictions with power corrections and in the Rgap scheme. 
The dashed vertical lines indicate the PDG 2014 [231 deter¬ 
mination of a s (mz)- 


that of Ref. (10 (higher statistics for the two-loop non¬ 
singular cross sections and using the exact result for the 
two-loop soft function non-logarithmic constant). In ad¬ 
dition we have corrected the systematic uncertainty for 
the ALEPH data, Q = 91.2 GeV of Ref. [64]Q When 
we compare thrust and C-parameter we neglect bottom- 
mass and QED effects in both event shapes. In this setup, 
we find an updated result for thrust: 

a s {m z ) = 0.1134 ± 0.0002 exp (35) 

± 0.0005hadr±0.0011pert, 

^i(Ra,Ha) = 0.329 ± 0.009 eX p 

± 0.021 Qa ( mz ) ± 0.060 P ert GeV. 

For completeness we also quote an updated thrust result 
when both QED and bottom-mass effects are taken into 
account: 

a s [m z ) = 0.1128 ±0.0002 exp (36) 


5 In Ref. 0 we assumed that two quoted uncertainties where asym¬ 
metric uncertainties, but it turns out they are two sources of sys¬ 
tematic uncertainties that need to be added in quadrature. This 
has no significant effect on the results of Ref. [9]. 


FIG. 14. Distribution of best-fit points in the a s (m z )- 2fR 
plane for both thrust (blue) and C-parameter (red) at N 3 LL'-|- 
0(a^)+Qi(R, fi). The outer solid ellipses show the Ay 2 = 2.3 
variations, representing l-tr uncertainties for two variables. 
The inner dashed ellipses correspond to the l-er theory uncer¬ 
tainties for each one of the fit parameters. The dotted ellipses 
correspond to Ay 2 = 1 variations of the total uncertainties. 
All fits have been performed with N 3 LL' theoretical predic¬ 
tions with power corrections and in the Rgap scheme. This 
plot zooms in on the bottom two ellipses of Fig. |13| 



a's(mz) 

FIG. 15. Distribution of best-fit points in the a s (m z )- 2fli 
plane for both thrust (blue) and C-parameter (red) at N 3 LL'-|- 
G(a 3 ) + fir. The meaning of the different ellipses is the same 
as in Fig. [14] 
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± 0.0005hadr±0.0011p ert , 
^i(Ra^a) = 0.322 ± 0.009exp 

± 0.021 as ( mz ) ± 0.064 pert GeV. 


a short-distance Rgap scheme, in which the sensitivity to 
infrared physics is reduced. 

As our final result from the C-parameter global fit we 
obtain, 


Both the results in Eqs. ([35]) and (36) are fully compatible 
at l-cr with those in Ref. |9j, as discussed in more detail 
in App. [Cj 

When testing for the universality of IR between thrust 
and C-parameter, there is an important calculable nu¬ 
merical factor of 37 t/2 = 4.7 between fij and ilf that 
must be accounted for; see Eq. pt). If we instead make 
a direct comparison of f l\ and Ilf', as shown in Fig. |l3 


(lowest blue ellipse vs uppermost green ellipse, respec¬ 
tively) then the results are 4.5-cr away from each other. 
Accounting for the 37r/2 factor to convert from Ilf to 
the upper green ellipse becomes the centermost red 
ellipse, and the thrust and C-parameter determinations 
agree with one another within uncertainties. Due to our 
high-precision control of perturbative effects, the fR pa¬ 
rameters have only ~ 15% uncertainty, yielding a test of 
this universality at a higher level of precision than what 
has been previously achieved. 

A zoomed-in version of this universality plot is shown 
in Fig. [14] The upper red ellipse again shows the re¬ 
sult from fits to the C-parameter distribution, while the 
lower blue ellipse shows the result from thrust tail fits. 
For both we show the theory uncertainty (dashed lines) 
and combined theoretical and experimental (dotted lines) 
39% CL uncertainty ellipses, as well as the solid ellipses 
which correspond to Ay 2 = 2.3 which is the standard 
l-cr uncertainty for a two-parameter fit (68% CL). We 
see that the two analyses are completely compatible at 
the l-cr level. An important ingredient to improve the 
overlap is the fact that we define the power corrections 
in the renormalon-free Rgap scheme. This is shown by 
contrasting the Rgap result in Fig. [14] with the overlap 
obtained when using the MS scheme for Hi, as shown in 
Fig. [15} 


VIII. CONCLUSIONS AND COMPARISON TO 
OTHER q s DETERMINATIONS 

In this paper an accurate determination of a s from fits 
to the C-parameter distribution in the resummation re¬ 
gion was presented. We fit to the tail of the distribution 
defined by 37 tAqcd/Q <C C < 3/4, where the dominant 
hadronization effects are encoded in the first moment of 
the shape function fR, which is a power correction to 
the cross section. By fitting to data at multiple Q’ s, 
the strong coupling a s (mz ) and fR can be simultane¬ 
ously determined. The key points to our precise theo¬ 
retical prediction are: a) higher-order resummation ac¬ 
curacy (N 3 LL'), achieved through an SCET factorization 
theorem, b) 0(a 3 ) matrix elements and fixed-order kine¬ 
matic power corrections, c) field-theoretical treatment of 
nonperturbative power corrections, and d) switching to 


a s (m z ) = 0.1123 ±0.0015, 
ni(R A ,PA)= 0.421 ± 0.063 GeV, 


(37) 


where a s is defined in the MS scheme, and Hi in the 
Rgap scheme (without hadron-mass effects) at the refer¬ 
ence scales = Ha = 2 GeV. Here the respective total 
l-cr uncertainties are shown. The results with individ¬ 
ual l-cr uncertainties quoted separately for the different 
sources of uncertainties are given in Eq. (34). Neglect¬ 


ing the nonperturbative effects incorporated by Hi, the 
fit yields a s (mz) = 0.1219 which exceeds the result in 
Eq. (37) by 8%. This is consistent with a simple scaling 


argument one can derive from experimental data, pre¬ 
sented in Ref. m We have also presented an updated 
thrust result, using our improved profiles for thrust and 
including bottom-mass and QED effects (but neglecting 
hadron-mass effects). This global fit for thrust gives 

a s {rn z ) = 0.1128 ± 0.0012, (38) 

Hi(J?A,A* a) = 0.322 ± 0.068GeV. 


Our theoretical prediction is the most complete treat¬ 
ment of C-parameter at this time, and, to the best of 
our knowledge, all sources of uncertainties have been in¬ 
cluded in our final uncertainty. Possible improvements 
which are expected to be negligible relative to our final 
uncertainty include finite bottom-mass effects, QED ef¬ 
fects, and axial-singlet contributions. From our results 
there are a number of theoretical avenues that lead to 
small effects but which would be interesting to investigate 
further in the future. These are common to almost ev¬ 
ery event-shape analysis in the literature and include (i) 
resummation of logarithms for the nonsingular partonic 
cross section; (ii) the structure of nonperturbative power 
corrections for the nonsingular contributions (the last two 
points can be clarified with subleading SCET factoriza¬ 
tion theorems); (iii) analytic perturbative computations 
of the G(aig) nonlogarithmic coefficients in the partonic 
soft function and the jet function, as well as the four-loop 
QCD cusp anomalous dimension (and to a lesser extent, 
the numerically determined constant of the two-loop 
partonic soft function); (iv) a better understanding of 
hadron-mass effects, and in particular their resummation 
beyond NLL; (v) a better theoretical description of the 
region around and above the shoulder. Concerning (i), 
and following the common lore, we have incorporated in 
our analysis the nonsingular contributions in fixed-order 
perturbation theory. However we have estimated the un¬ 
certainty related to the higher-order logarithms through 
the usual renormalization scale variation. Concerning (ii) 
we observe that the effect of these neglected power cor¬ 
rections is much smaller than naively expected, as can 
be seen from a comparison of our theoretical prediction 
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and LEP data in the far-tail region. A first step towards 
clarifying (i) and (ii) has been taken in Refs. [ZH[75], for 
the case of thrust. The computation of missing pertur¬ 
bative terms (iii) is a priori feasible with current com¬ 
putational knowledge but they do not dominate our per¬ 
turbative uncertainties. Concerning (iv) we have shown 
that hadron-mass effects have a very small impact on 
the determination of a s , and hence unless the rest of the 
sources of uncertainty become substantially smaller, our 
lack of knowledge does not constitute a problem. As for 
(v), our fits do not include data above the shoulder, so 
this problem has no impact on our fit. Nevertheless an 
analysis of these subleading effects would be interesting. 

The same theoretical program carried out for thrust 
and C-parameter can be applied to other event-shapes, 
and the most prominent one is Heavy-Jet-Mass. This 
has partially worked out already in Ref. [8] at the purely 
perturbative level using fully canonical profiles. Their de¬ 
termination of a s is discussed below. For recoil-sensitive 
observables such as Jet Broadening [75H5D] , one needs to 
deal with rapidity singularities, which imply that addi¬ 
tional logs need to be resummed, and more complicated 
nonperturbative power corrections. The former has been 
pushed to the N 2 LL order in Ref. [51] . and the latter has 
been studied in Ref. [[82]. Recoil-insensitive versions of 
Broadening have also been derived [S3], but not yet stud¬ 
ied experimentally. Finally, it is very straightforward to 
generalize our theoretical treatment to the case of ori¬ 
ented event shapes [84] . in which one additionally mea¬ 
sures the angle between the beam and the thrust axes. 

At this point we compare our result for a s with other 
determinations from event shapes at O(a^). To the best 
of our knowledge, the only analyses which fit to the tail 
of the C-parameter distribution using three-loop input 
are Ref. [55] (using purely fix-order perturbation the¬ 
ory) and Ref. [86] (including NLL resummation). Both 
analyses use Monte Carlo (MC) event generators to es¬ 
timate hadronization effects, and fit a s for different Q 
values, finding values a s (mz) = 0.1288 ± 0.0043 and 
0.1252 ± 0.0053 respectively for a fit to the Q = 91.2 GeV 
data. These larger a s (mz ) values are consistent with 
our fits which neglect power corrections, and follow¬ 
ing Ref. [1] we can conclude from this that MCs does 
not provide a reasonable estimate of the power cor¬ 
rections when including the higher-order perturbative 
contributions. In Ref. [55] two-parameter global fits 
to the first five moments of the C-parameter distribu¬ 
tion were performed. Hadronization effects are included 
via the frozen coupling model, and the value obtained, 
a s (rnz) = 0.1181 ± 0.0048, is fully consistent with our 
result in Eq. (371 at 1 -a. 


A graphical comparison with other event-shape deter¬ 
minations is shown in Fig. |16| The figure includes de¬ 
terminations where power corrections are estimated with 
MC generators, labeled by 7-10. Analyses 1-6 correspond 
to those in which power corrections were incorporated 
with an analytic method (either a shape function or the 
dispersive model). In the analyses 1-6 global fits are per- 


n —i— 1 — 1 — 1 — 1 —i— 1 — 1 — 1 — 1 —i— 1 — 1 — 1 — 1 —i— 1 — 1 — 1 —■—i— 1 —■— 1 — 1 —i—■— 1 — 1 —■—i— 1 —■—■■ 

Analytic power corrections, moments + 0(a)) 

1 _ Thrust Moments 

(AFHMS '12) 

2 ,_ @ _, t, p, Y 3 , B Moments 

(GJL '09) 

Analytic power corrections, distribution + Ota)) 

3 ,_•_, C-parameter N 3 LL 

(HKMS '14) 

4 Thrust N 3 LL 

(AFHMS '10) 

5 m Thrust n2ll 

(GLM’12) 

Thrust NLL 

6 ' * 1 (DW '09) 


MC power corrections, distribution + 0(a)) 

7 HJM N 3 LL 

(CS '10) 

8 Thrust N 3 LL 

(BS ’08) ' * ' 

9 r,p,Y 3 ,B,NLL ,_._ 

(DGGGHLS '09) 

10 r, p, Y 3 , B w/o resum. ,_,_, 

(DGGGHS '08) 

_ 1 ■■■ 1 ... 1 ... 1 ... 1 ... 1 1 • 

0.110 0.115 0.120 0.125 

a s (m z ) 


FIG. 16. Comparison of our determination of a s (mz ) (red) 
with similar analyses from thrust [2 [TO] (green) and other 
determinations from fits to event-shape distributions using 
(9(a)) theory predictions and different levels of resummation. 
Results shown below the lower dashed line include power cor¬ 
rections as predicted by MC generators, and results above 
this line treat power corrections either from a shape func¬ 
tion (red and green) or from the dispersive model (orange Q] 
and purple ED- Determinations above the upper dashed line 
correspond to fits to moments of the distributions, and those 
below to fits to the tail of the differential distribution. The 
translucent green band corresponds to the world average from 
Ref. [HH- 


formed, whereas in the 7-10 analyses a s was determined 
at multiple Q values and the final result is an average of 
those. Only analyses 1 [TO], 3 (this work), and 4 [5] used 
a completely field-theoretical approach for the power cor¬ 
rections. We also show both results from fits to the event- 
shape distributions (3-10) and from fits to moments of the 
event shape distribution (1 and 2). Although all analy¬ 
ses included 0 (a )) matrix elements, different levels of 
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tors; a determination analyzing the lattice prediction for 
the QCD static potential |29j; a reanalysis of electroweak 
precision observables by the Gfitter collaboration E3; 
the most recent analysis of tau decays in which the re¬ 
cently released ALEPH data was used together with the 
OPAL data; the previous determinations from fits to the 
thrust distribution [9] and moments of the thrust distri¬ 
bution da; and of course the current world average [23] 
(shown as the green band). The ABM (DIS) and thrust 
results are compatible with our determination, while in 
contrast the disagreement with either lattice QCD or the 
world average is 4 -a. Many other determinations lie be¬ 
tween these two values. The source of this disagreement 
is an important open question. 


FIG. 17. Comparison of our determination of a a (mz) 
(red) with similar previous analyses at N 3 LL' for thrust 
(green) [9J [TU] and other selected determinations: lattice [24] 
and static energy potential [29] (both use lattice input, in 
blue), Electroweak precision observables fits [87] (black), Deep 
Inelastic Scattering [25] and global PDF fits [501 I5T1 . and 
hadron r decays [27] (Fixed Order Perturbation Theory lower, 
and Contour Improved Perturbation Theory, both in gray). 
The current world average [23] is shown as a translucent green 
band. 


resummation have been achieved. Analyses 2 [25] and 
10 [55] did not included resummation; 6 [4] and 9 [SB] 
included NLL resummation; 5 m include N 2 LL resum¬ 
mation; and analyses 3,4,7, and 8 included N 3 LL resum¬ 
mation. Analyses 2 , 9, and 10 simultaneously fit to many 
event shapes, whereas the others focused on a single ob¬ 
servable: thrust (1, 4-6 and 8 [3]), Heavy-Jet-Mass (7 [5]), 
and C-parameter (3, which is this work). The analyses 
1, 3, 4, 7 and 8 used SCET to perform Sudakov log re¬ 
summation. All results that used an analytic treatment 
of power corrections have smaller values of a s . This is 
consistent with a simple dimensional analysis argument 
(see Refs. ® H 2 )- Higher order resummation results 
in a convergent perturbation series and smaller uncer¬ 
tainties, and the Rgap scheme also reduces uncertain¬ 
ties. Accounting for the fact that results relying on MC 
for the treatment of power corrections should likely have 
larger hadronization uncertainties, all results are com¬ 
patible among one another. The most precise results are 
however clearly in disagreement with the world average, 
which is dominated by lattice QCD results (see below) 
and shown as a translucent green band. 


We conclude this work by comparing our result for 
ct s (rnz ) with the results of a selection of recent analy¬ 
ses usin g other techniques and observables, as shown in 
Fig. 17 We include a N 3 LO analysis of data from deep 
inelastic scattering from the ABM group [28], the global 
PDF fits of the MSTW group [31] and the NNPDF col¬ 
laboration [30]; the most recent (and accurate) determi¬ 
nation from the HPQCD lattice collaboration [22j, from 
the analysis of Wilson loops and pseudoscalar correla- 
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Appendix A: Profile Formulae 


In this appendix, we give the details for the profile 
functions that control the renormalization scales as laid 


out in Sec. IIB For the soft profile function, we use the 
form, 


Ms = 


Mo 

0 < C < t 0 

C(Mo, 0, 0, to, h, C) 

t 0 <C < t\ 

r s Vh § 

t\ < C <c £2 

C( 0 , /v/,o, t2,ts,c) 

t 2 < c < t s 

M H 

t s < c < 1 


(Al) 


where the physical meaning of the parameters is ex¬ 
plained in Sec. II B| The function £(ai, &i, 02, 62, ti, £2, t) 
(with t\ < £ 2 ), which smoothly connects two straight 
lines of the form h(t) = aq + b\ t for t < t\ and 
12 (f) = cl 2 + 62 f for t > t 2 is given by 


C(*) = 


fii + bi(t - ti) + ei (t — <i) 2 

ci 2 + b 2 (t - t 2 ) + e 2 (i - ^2) 2 


tl — t — tm 
tm — t — ^2 
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FIG. 18. Comparison of a s determinations from C-parameter 
tail fits in the thrust Rgap scheme (lower red ellipse) and the 
C-parameter Rgap scheme (upper blue ellipse). The leading 
power correction fi) in the C-parameter Rgap scheme is con¬ 
verted to fii in the thrust Rgap scheme in order to have a 
meaningful comparison. Theoretical uncertainty ellipses are 
shown which are suitable for projection onto one dimension to 
obtain the 1-cr uncertainty, without experimental uncertain¬ 
ties. 


di — ai + b\ t\ , q >2 — (i 2 + &2 t 2 5 

_ 4 (5,2 — ai) — (3 b\ + 62 ) (^2 ^ ti) 

61 “ 2(f 2 -G) 2 

_ 4 (ai — a 2 ) + (3 b 2 + 61 ) (t-2 — ti) 

62 “ 2 {t-2 - ti) 2 

For the jet scale, we use the form 


(A2) 


MC9 


[1 + e J {C-t a ) 2 ] s /n H Hs{C) C<t s 
M h C > t s 


(A3) 


which allows a slight modification of the natural relation 
between the scales /ij = i n order to account for 

theoretical uncertainties. 

For the subtraction scale, we have 


) Ro 0 < C <C to 

C(f?o, 0, 0 , to, h,C) t 0 < C < ti ■ 
fis(C) h < C < 1 

( A4 ) 

As explained earlier, we take R = /is in the resum¬ 
mation region to avoid large logs and R ^ /is in the 
nonperturbative region to remove the renormalon. The 
C, function here interpolates smoothly between these two 
regions. 


FIG. 19. Comparison of thrust a s determinations using our 
new profiles (left red ellipse) and the profiles of Ref. [9j (right 
blue ellipse). Theoretical uncertainty ellipses are shown which 
are suitable for projection onto one dimension to obtain the 
1 -cr uncertainty, without experimental uncertainties. 


It is necessary to vary the profile parameters to es¬ 
timate the theory uncertainty. We hold the difference 
between the parameters associated with the purely non¬ 
perturbative region constant: /r 0 — Rq = 0.4 GeV, and 
we set as default values ^0 = 1-1 GeV, Ro = 0.7 GeV. 
We are then left with nine profile parameters to vary 
during the theory scan, whose central values and vari¬ 
ation ranges used in our analysis are: r s = 2 x 1.13 ±:l , 
n 0 = 12 ± 2, m = 25 ± 3, t 2 = 0.67 ± 0.03, 
t s = 0.83 ± 0.03, ej = 0 ± 0.5, en = 2 ±:l and n s = 0±1. 
These variations are shown in Tab. |U 


Appendix B: Comparison of thrust and C-parameter 
subtractions 


In Fig. [18] we compare fits performed in the Rgap 
scheme with C-parameter gap subtractions as the upper 
red ellipse, and for our default fits in the Rgap scheme 
with thrust gap subtractions as the lower blue ellipse. 
At N 3 LL/ order with C-parameter subtractions the re¬ 
sults are a s (mz) = 0.1126 ± 0.0002 exp ± 0.0007hadr ± 
0.0022 pert and fii(i?A,A*A) = 0.447±0.007 exp ±0.018« s ± 
0.065 pert GeV, with x 2 lin /dof = 0.988. One can see that, 
even though both extractions are fully compatible, the 
thrust subtractions lead to smaller perturbative uncer¬ 
tainties. This is consistent with the better perturbative 
behavior observed for the cross section with thrust sub¬ 
tractions in Ref. (T2] . 
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Appendix C: Comparison of thrust results with 
Ref. [9] 

In Fig. [19] we compare global fits for the thrust distri¬ 
bution using the profiles of Ref. [9] (shown by the right 
ellipse in blue) and the profiles used here (shown by the 
left ellipse in red). As mentioned earlier, the profiles 
used here have several advantages over those of Ref. [2] 


in terms of their ability to independently impact the dif¬ 
ferent regions of the thrust distribution, and in particular 
do a better job in the nonperturbative region (which is 
outside our fit region). The two versions of the profiles 
are consistent within their variations, and the fit results 
shown for 39% CL for two dimensions in Fig. [l9] (which 
is 68% CL for each one-dimensional projectionjare fully 
compatible. 
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